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ABSTRACT 

We study how joint shear and magnification measurements improve the statistical precision of weak lensing 
mass calibration experiments, relative to standard shear-only analysis. For our magnification measurements, 
we consider not only the impact of lensing on the source density of galaxies, but also the apparent increase in 
source sizes. The combination of all three lensing probes — density, size, and shear — can improve the statis- 
tical precision of mass estimates by as much as 40% -50%. This number is insensitive to survey assumptions, 
though it depends on the degree of knowledge of the parameters controlling the magnification measurements, 
and on the value of the parameter q that characterizes the response of the source population to magnifica- 
tion. Furthermore, the combination of magnification and shear allows for powerful cross-checks on residual 
systematics (such as point spread function corrections) at the few percent level. 

Subject headings: cosmology: clusters - weak lensing 
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1. INTRODUCTION 

Estimating cluster masses is a key component of many 
cosmological and astrophysical experiments. Cosmologi- 
cally, the abundance of galaxy clusters as a function of their 
mass and redshift is widely recognized as a powerful probe 
of the fundamental p hysics behind our accelerating universe 
dAlbrecht et alj|2006l) , and currently provides one of the most 
robust estimates of the amp litude of the low r edshift mat- 
ter power spectrum trg (e.g. iMantz et al J 120081 iHenrv et"afl 
l2009HVikhlinin et al.ll2009t lRozo et alj|2010al) . Clusters have 
also been used to con s train gravity on cosmological scales 
dSchmidt et al.l I2009bt iRapetti et all 120 10t lLombriser etafl 
12.010b . Astrophysically, comparison of weak lensing mass 
estimates to those derived from X-ray or SZ observables 
can be used to study halo structure and the importance 
of non-thermal processes in the intra-cluster medium (e.g . 



Mahdavi et alJ2008l:)Marrone et al.l2009tlUmetsu et all2009t 



Morandi et all 120 lot iMeneghetti et al.l 12010b iMolnar et at] 
20101) . Moreover, precise cluster mass estimation signifi- 
cantly enhances the utility of galaxy clusters as a probe of 
galaxy evolution by allowing us to select galaxy populations 
in halos of known masses. The fact that all these applica- 
tions require robust and well understood cluster mass esti- 
mates highlights the importance of cluster mass estimation as 
a fundamental tool of contemporary astronomy. 

Today, weak gravitational lensing is widely identified as 
one of the most promising t ools for robustly estimating clus - 
ter masses (for a review, see lBartelmann & Schneiderll2001l) . 
In the standard approach, one relies on the fact that grav- 
itational lensing by foreground structures will shear back- 
ground galaxies, introducing a tangential pattern in the ellip- 
ticities of background galaxies about the center of the fore- 
ground overdensity. These, however, are not the only ob- 
servable signatures of weak gravitational lensing, as the lat- 
ter also impacts the number density, magnitudes, and sizes 
of background galaxies. In this letter, we explore how the 
combination of all three lensing measurements can help esti- 
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mate cluster masses more accurately and robustly by break- 
ing degeneracies between the desired physical quantities and 
systematic nuisance parameters. While galaxy s izes has al - 
ready been proposed as a weak lensing signal (e.g. lJainl2002l) . 
and joint shear and density analysis has already been ex- 
plored in the litera ture — both in the context of cosmic shear 
dZhang et alJ2010h a nd for weak lensing mass calibration (for 
some examples see iBroadhurst et ail 119951: ISchneider et alj 
120051 IVanWaerbeke et all 120 lot lUmetsu et al.l 1201 Oh — we 
consider for the first time all these observational signatures 
simulteneously. Moreover, we do not ignore the covariance 
between the various observational signatures of gravitational 
lensing, as has been done in previous studies. 

As we were finishing this letter, we learned of an indepen- 
dent analysis that is very much on the same spirit as ours, but 
which focuses on cosmic shear, and in particular on the ability 
of jo int lensing probes to se lf-calibrate multiplicative shear er- 
rors (IVallinotto et al.l2010T) . The two analysis are qualitatively 
similar and reach very similar conclusions, namely that joint 
lensing analysis are a promising avenue of self-calibrating 
systematic uncertainties in weak lensing experiments. 

2. THE WEAK LENSING OBSERVABLES 

We consider the problem of weak lensing mass calibration 
using three different observables. All of our estimators use ra- 
di al bins r n , and globa l normalization. We adopt the notation 
of iRozo et alJ d2010bl) . The first observable is n a , the mean 
source density in an annulus: 



(1) 



In the above expression, we have pixelized all space, the index 
i runs over all spatial pixels, «, is the source density at pixel i 
(so that «, AS1 is either or 1), W" is the top hat filter defining 
the annulus of interest, and A a is the area of said annulus. 

The second observable we consider is the area-averaged 
number weighted shear of galaxies in an annulus, which we 
denote h. Specifically, we write 



A a 



where e, is the galaxy ellipticity field at pixel i. 



(2) 
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Finally, the third observable is the number-weighted aver- 
age apparent area of source galaxies s a , defined via 



1 

A a 



(3) 



where a -, is the area of the source galaxy under consideration. 
Note that one can choose among different po ssible estima- 
tors to capture magnification via galaxy sizes 
Generally, the different estimators have similar statistical and 
systematic properties, and our conclusions do not depend on 
the precise choice of estimator. 
For simplicity, we assume sources are unclustered so that 
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where q is the parameter governing the impact of lensing bias, 
i.e. n bs = pfil 2 n tme dSchmidt et al.ll2009ah . and fa is the mag- 
nification in the direction of galaxy ;. Source clustering will 
be briefly discussed in section|5] We further assume the ellip- 
ticities of the sources satisfy 



?i e j)=gigj + <>ij cr e/ 2 - 



(5) 
(6) 



As for the source area, it is obviously of critical importance 
to consider the impact of the Point Spread Function (PSF) on 
the apparent galaxy area. For Gaussian sources and PSF, the 
PSF adds to the apparent size, so we assume 



@obs ~ fl&true ^psf 



(7) 



where «p S f characterizes PSF broadening and where we have 

included the impact of magnification on the apparent area 
of the source. Usually, size estimates are corrected (decon- 
volved) for PSF effects, so in principle, the PSF offset a^ s f 

is identically zero. In practice, imperfect PSF deconvolution 
can lead to systematics offsets, so we will retain a^ s f as a 

systematics parameter. 
Defining 

a = (a true ) (8) 
a fl 2 = Var(<w)/a 2 (9) 
e = a psf/ 5 ( 10 ) 

(ai)=a(p + e) (11) 
(a i a j ) = {a i ){a j ) + 5 i jt J L 2 a 2 a 2 a (12) 

with e characterizing a possible PSF correction bias. 

With these assumptions, the expectation values and corre- 
lation matrix of our observables are given by 

1 -nut 



one has that 
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where a,b £ {n,h,s}, and the symbol ~ signifies that the co- 
variance matrix is symmetric. In deriving the above equations 
we have made use of the narrow-bin approximation, that is, 
we have assumed all lensing quantities are constant across the 
radial bins employed for the measurements. 

It is useful to allow for the possibility of the source density 
of galaxies used for the density measurement n to be different 
from the source density used for the shear and area measure- 
ments. In such a scenario, the source galaxies employed for 
the density measurement that are not included in the source 
density measurement will in general have a different value for 
the q parameter. Assuming the relative densities are related by 
a factor A such that the source density for measuring n is A«, 
and letting Q denote the lensing parameter of the sources em- 
ployed in the density measurement, the above equations are 
modified only via 
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(17) 
(18) 

(19) 
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Note that the covariance between density and other observ- 
ables goes as yfl' 2 since only galaxies that are shared by the 
two measurements are employed, so the appropriate lensing 
parameter is q rather than Q. 

We will be primarily interested on how a joint analysis of 
the above lensing observables compares to the standard weak 
lensing shear analysis, which uses as its observable the mean 
galaxy ellipticity in annuli, 



E,.Ar>n, e ,W, a 



(21) 



Note that this estimator uses a location-dependent normaliza- 
tion, rather than a global norm alization. The mean and vari- 
ance of this observable are (see lRozo et aT] ( l2010bl) ) 
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We use the above set of equations to set the baseline against 
which the results from a joint analysis are to be compared 
against. 

2.1. Signal-to-Noise Considerations 

Before we forecast the precision with which weak lensing 
observables can constrain cluster masses, it is worth consider- 
ing the signal-to-noise of each of our three observables: den- 
sity («), shear (h), and size (s). For simplicity, we will work 
in the weak lensing limit and assume negligible systematics 
(e = 0). We find 



(S/N)„ = (nA)^ 2 X^ 2 QK 
(S/N) h = (iiA)^ 2 V2^ 

(S/N) s = (nA) l/: 



[1 + * 2 ]V 2 



(24) 
(25) 

(26) 
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There are two key points to take away from the above re- 
sults: 

1. The signal-to-noise of all measurements scales as 
(nA) 1 / 2 . Consequently, the ratio of the mass uncer- 
tainty between any two experiments is independent of 
the assumed source density. We use this to our advan- 
tage by normalizing all of our uncertainties relative to 
that of a standard shear-only experiments. Note, how- 
ever, that when including priors our quantitative results 
are not entirely «-independent since the relative impor- 
tance of priors varies with n. That said, we find this n 
dependence on the relative errors to be rather modest 
for moderate changes in h. 

2. The signal-to-noise of the size measurement is indepen- 
dent of a. Consequently, all of our results are fully in- 
dependent of the fiducial value of a, a fact which we 
have explicitly checked. 

One last important result that can be garnered from the 
above signal-to-noise estimates is that large boost parameters 
A are not necessarily optimal for density measurements. For 
instance, in current experiments, if source selection occurs 
close to the confusion limit then A can be as large as A ps 4, 
but the Q parameter tends to be relatively small, Q « 0.5. On 
the other hand, if one requires that only sources well above 
the confusion limit be included in the source population, as 
is done for shear and area measurements, then A is modest, 
with 1 < A < 1-5, but Q can be significantly larger, with 
Q « q w 1.5 (ISchmidt et al .11200 9 a). Comparing the former 
choice with (Q, A) = (0.5,4) to the latter case of (1.5, 1.5), we 
find that the second choice is expected to have nearly twice the 
signal-to-noise. Consequently, from now on we will make the 
simplifying approximation that Q = q and the boost factor A is 
moderate, no more than 1.5. In practice, given an empirically 
determined relation Q{X), one may find the source cuts which 
maximize the signal-to-noise of the density measurement by 
setting 

^5=0. ,27, 

In principle, if Q keeps increasing quickly with decreasing 
A, it is possible that the density measurement is optimized 
using fewer source galaxies than those employed in the shear 
measurement. 

3. FIDUCIAL MODEL AND FISHER MATRIX FORECAST 

We assume halo m ass pro files are described by the standard 
iNavarro et al.l (NFW. [l996l) form, which is determined by two 
parameters, the halo mass M200m> defined relative to the mean 
matter density, and its concentration. We set M200/11 = 10 15 M Q 
and c = 5, and adopt lens and source redshifts of zl = 0.3 and 
z s = 1-0 respectively. For our background cosmology we as- 
sume flat ACDM model with Sl,„ = 0.28 and h = 0.7. Distances 
are in physical units in Mpc, not /r'Mpc. 

We assume each of our three observables — n, h, and s — 
is measured in narrow, logarithmically spaced radial bins of 
width A \nR = 0.05, extending from R m j„ = 0.2 Mpc to R„ mx = 
2 Mpc. Our results are largely insensitive to A In/?, and only 
mildly sensitive to the innermost radius /?„„„ so long as the 
strong lensing region is excluded. There is some sensitivity 
to the assumed maximum radius R max , in the sense that larger 
Rmax values decreases the relative importance of density and 
size measurements. 



TABLE 1 
Fisher Matrix Parameters 



Parameter 


Fiducial Value 


tt 1 

Explanation 


M 


10 15 M Q 


Halo mass 


c 


5 


Halo concentration 


n 


10 sources/arcmin 2 


Mean source density 


<i 


1.5 


h and s lensing parameter 


a 


Irrelevant 


Mean source size 


e 


0.0 


Systematic Size PSF Residual 


<Je 


0.4 


Std. Dev. of source ellipticity 


da 


1.2 


(Std. Dev. of source area)/a 


A 


1.5 


Density measurement source boost 


ZL 


0.3 


Halo redshift 


Zs 


1.0 


Source redshift 


Note. — 


Quantities below the horizontal lines are held fixed and 



assumed to be known a priori throughout the entire manuscript. 

For the lensing parameters, we set q = 1 .5 and e = 0. We 
also adopt a e = 0.4 and a a = 1.2 for the intrinsic ellipticity 
and size dispersion. The former value is typical of ground- 
based shear experiments, while the latter is estimated from an 
ongoin g study of size-lensin g with data from the COSMOS 
survey dScoville et al.ll20O7l) . In accordance with the discus- 
sion in section |2~T1 we force the parameter Q to be identical 
to q — i.e. we assume the same size and magnitude cuts in 
source selection for all measurements — and we set A = 1.5. 
For the source density, we assume h = 10 galaxies/arcmin 2 , 
as expected for the DES, though we note our results are fairly 
insensitive to the precise value of n. The fiducial value for a 
is irrelevant as our results are a-independent (see section lXTI ). 
That said, the precision with which a is known is important, 
so we do consider how uncertainties in a impact our results. 

The Fisher matrix for a joint density+shear+size weak lens- 
ing experiment is given by 



d(a a ) 



a a,b 



dpi dpj 



(28) 



where a,b loop over our three observables, and a indexes the 
radial bins employed in the experiment. The vector p is the 
vector of parameters of interest, which at minimum includes 
halo mass and concentration, but can also include additional 
nuisance parameters such as the mean source density n or the 
lensing bias parameter q. The corresponding covariance ma- 
trix is given by C = F~ l . Table Q] summarizes the parameters 
considered in our Fisher Matrix analysis. 

4. RESULTS 

Figure [T] shows our forecasted 68% error ellipse in the 
mass-concentration plane for a variety of weak lensing mass 
calibration experiments. The mass and concentration axes 
have been displaced so that the input parameters occur at the 
origin, and the axis are in units of the marginalized standard 
deviation of each parameter for a standard shear analysis. 
When normalized in this fashion, our results are insensitive 
to the source density n. All parameters ot her t han halo mass 
and concentration are held fixed (section 14.11 will relax this 
assumption). 

The top panel in Figure[T]considers each of our observables 
as an independent weak lensing experiment. The solid line 
corresponds to the standard shear-only analysis, and is nearly 
identical to that obtained using the observable h alone (not 
shown). The dotted and dashed lines correspond to density («) 
and size (s) respectively. We find shear measurements lead to 
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FIG. 1. — 68% error ellipses in the halo mass-concentration plane for a 
variety of weak lensing experiments as labelled. The top panel treats each 
observable independently while the bottom panel considers a joint analysis. 
The axes are chosen such that the input halo mass and concentration occur at 
the origin, and the scale of the axis is in units of the marginalized lcr error 
for a standard shear-only analysis. All parameters other than halo mass and 
concentration are held fixed. Joint analyses have the potential to increase the 
precision of weak lensing mass calibration experiments by 35%. 

the tightest error ellipses, but size and density measurements 
can achieve comparable and even higher precision than the 
shear measurement with respect to mass estimation. This con- 
clusion does depend on our fiducial assumptions, and a lower 
q value can render shear the most precise estimator. The im- 
portant point to emphasize though is that broadly speaking the 
various estimators have comparable precision. 

The bottom panel in Figure Q] illustrates how the pre- 
cision of a shear experiment improves as we include ad- 
ditional observables. The curves shown are for a shear 
only (solid), shear+density (dotted), shear+size (dashed), and 
shear+size+density experiment (dash-triple dot). We find that 
a joint analysis of all three observables significantly improves 
the precision of weak lensing mass calibration experiments. 
Specifically, the uncertainty in the best fit log-mass is reduced 
by almost 50%, with cr lnM j i nt / 

°inM shear only — 0.53. From 
now on, whenever we quote mass uncertainties, they will al- 
ways be normalized in this way, so that <t\ b m = 1 mean that the 
analysis under consideration achieves the same precision as a 
shear-only analysis. 

4.1. Priors and Systematics Self-Calibration 



FIG. 2. — 68% confidence contours for a joint shear+density+size weak 
lensing mass calibration experiment. The origin corresponds to the input halo 
mass and concentration, and the scale of the axes is in units of the marginal- 
ized lcr eiTor for a standard shear-only analysis. The solid ellipse is the stan- 
dard shear-only results, while the remaining ellipses are for a joint analysis 
with various priors. They are as follows. Dashed: fixed n, q, and a, with 
no prior on e, resulting in o\ n M = 0.59. Dotted: no priors on n, q, a, or e, 
resulting in cr{ n M = 0.93. Dash-triple-dot: fixed n, q, a, and e, resulting in 
o"i„m = 0.53. 

We have seen that joint shear+density+size analysis are 
in principle significantly more powerful than standard shear- 
only analysis. In practice, however, the various lensing pa- 
rameters that impact the density and shear signals are not per- 
fectly known a priori. We now explore to what extent our re- 
sults generalize when faced with uncertainties in n, q, a, and 
e. 

Figure |2] compares the 68% confidence contours of a stan- 
dard shear-only analysis (solid) to those of a joint analysis 
with a variety of priors for the nuisance parameters: perfect 
knowledge (dash-triple-dot), fixed n, q, and a but no prior on 
e (dashed), and no priors on any lensing parameter (dotted). 
Even in this last case, we find some small improvement in the 
precision of a weak lensing mass calibration experiment, with 
oinM = 0.93, i.e. a 7% improvement over the shear-only case. 

We want to determine now the level of precision required on 
a priori knowledge of our nuisance parameters in order to sat- 
urate the perfect knowledge limit illustrated in Figure|2] Since 
the density measurement is affected by only two of our four 
nuisance parameters, we begin by focusing on a shear+density 
experiment. 

The top panel in Figure[3]shows the error of a shear+density 
analysis as a function of the priors A Inn and Aq. In all 
cases, the error is normalized to that of a standard shear-only 
experiment. In these units, the minimum error achieved by 
a joint analysis when both n and q are perfectly known is 
cinM = 0.62, while complete ignorance of these parameters 
results in ai n M = 0.94. A priori knowledge of fi and q at the 
level of o~\ n n = 0.01 and a q = 0.1 is nearly as effective as hav- 
ing perfect knowledge of said quantities. Such level of a priori 
knowledge should be easily achievable. 

The solid curves in the bottom panel of Figure [3] trace the 
contours of the error in log-mass for our fiducial cluster in a 
joint shear+density+size analysis. We have assumed realistic 
priors Aln« = 0.01 and Aq = Q.l as per the above discussion, 
and we used lna as a parameter rather than a to enforce pos- 
itivity. With our adopted priors, the error achieved when In a 
and e are perfectly known is a\ n M = 0.56. To saturate this limit, 
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FIG . 3 . — Top panel: Error of a joint shear+density experiment as a function 
of the priors A Inn and Aq. The errors have been normalized to those of a 
standard shear-only analysis. In these units, perfect knowledge of n and q 
results in <T{ r m = 0.62, while a complete lack of knowledge results in oi nM = 
0.94. The perfect knowledge limit is nearly saturated with A Inn = 0.01 and 
Aq = 0.1. Bottom panel: Same as above, but now as a function of the priors 
A lna and Ae, assuming additional priors Aln/i = 0.01 and Aq = 0.1 (solid 
curves). The limiting values are cri nM = 0.56 when a and e are fixed, and 
"InM - 0.66 when neither quantity is known a priori. The perfect knowledge 
limit is nearly saturated at Alnd = 0.01 and Ae = 0.01. The dashed curve 
shows the a-posteriori error in e as a function of A In a when no priors are 
placed on e, assuming a source density it = 10 gals/arcmin - . Larger source 
densities will lead to tighter constraints in e. 

both lna and e must be known at the w 0.01 level, with the er- 
ror modestly increasing in the 0.01-0.1 range. Once lna and 
e are known to 10% level precision, further degradation oc- 
curs slowly. No knowledge of d and e results in a\ n M = 0.66. 
We expect measuring a should be no more difficult than mea- 
suring n, and thus Alna = 0.01 should be easily achieved. 
Saturating the e bound, on the other hand, is likely to prove 
difficult if not impossible, particularly for ground-based imag- 
ing. Thus, in practice we expect to be able to nearly saturate 
the limit in which e is left free, but «, q, and a are perfectly 
known (dashed ellipse in Figure|2]i. In this limit, o\ n M = 0.59, 
a 40% improvement over the standard shear-only analysis. 

There is an additional benefit to operating in this limit. 
The dashed curve in the bottom panel of Figure [3] shows the 
forecasted error in e for a joint density+shear+size analysis 
as a function of the prior Alna when no prior is placed on 
e. We see that even with the relatively low source density 
h = 10 sources/arcmin 2 we have assumed, a joint analysis can 



FIG. 4. — Error in log-mass as a function of halo mass for a joint 
shear+density+size weak lensing experiment, normalized to the standard 
shear-only result. Different curves correspond to different priors of the nui- 
sance parameters n (source density), q (lensing bias parameter), a (mean 
source area), and e (residual PSF correction to source area). The three cases 
we consider are no priors on nuisance parameters (dotted), no priors on e 
and ti (dash-dot), no priors on e (dashed), and all nuisance parameters fixed 
(solid). We expect real-world experiments will fall somewhere between the 
dashed and dash-dot curve. 

recover an unknown bias e to high precision. Assuming lna is 
known at the 1%- 10% level, the corresponding constraint 
in e is <r e sa 0.04-0.1. Note, too, that unlike the relative 
uncertainty between a joint lensing experiment and a shear- 
only analysis, the absolute error on e is in fact independent 
on the assumed source density. In particular, stacked weak 
lensing analysis, in which the effective source density can be 
a hundred or even a thousand times larger, will have a corre- 
spondingly tighter constraint on e. This is an incredibly useful 
result: if e is believed to be known a priori based on our un- 
derstanding of PSF subtraction, the joint analysis proposed 
here allows us to explicitly check whether or not e is indeed 
consistent with zero at the expected level. This is a highly 
non-trivial systematics check on PSF subtraction, particularly 
when lna is known a priori. 

5. SUMMARY AND DISCUSSION 

Gravitational lensing changes the density on the sky, size, 
and ellipticity of background galaxies. We have explored the 
extent to which a joint weak lensing mass calibration experi- 
ment that combines all three signatures can improve upon the 
standard shear-only approach. To perform this analysis, we 
must introduce 4 additional nuisance parameters: the mean 
source density n, the lensing parameter q governing the source 
density boost due to gravitational lensing, the mean source 
area a, and a PSF subtraction residual e that allows for size 
estimates to be systematically biased. 

Figure |4] summarizes our results. The figure shows the ra- 
tio between the uncertainty in log-mass for a joint analysis 
relative to the standard shear-only result. This ratio is plot- 
ted as a function of halo mass, and all errors are marginalized 
over all other parameters. The fiducial value for the halo- 
concentration as a function of mass is taken to be 

c(M)=5 -°(^y A (29) 

in agreement with the results of iNeto et alj ( 120071) . The vari- 
ous different curves correspond to different sets of priors: no 
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priors on nuisance parameters (dotted), no priors on e or h 
(dash-dot), no priors on e (dashed), and all nuisance parame- 
ters fixed (solid). 

In the limit that no nuisance parameters are known a priori, 
a joint analysis provides only a marginal improvement relative 
to the shear-only case (dotted curve). In practice, however, 
we expect the lensing parameters «, q, and a to be known 
with sufficient precision to saturate the limit in which these 
quantities are perfectly known (dashed line). For such a sce- 
nario, the precision of weak lensing mass calibration experi- 
ments can improve by as much as 30% relative to the standard 
shear-only case. Further improvements are possible if e can 
be controlled at the few percent level, which is likely to be a 
challenging demand, at least for ground-based surveys. 

Throughout, we have not incorporated source clustering 
in our error estimates. While in principle source cluster- 
ing affects all lens ing measurements (see the discussion in 
iRozo et ail (1201 Obi) '), the density measurement is affected di- 
rectly and most significantly. Source clustering will be impor- 
tant if ^{6)(nA) l l 2 > 1, where £(f?) is the source galaxy angular 
correlation function, and 9 = \J rAr/ cIa(zl) is the mean source 
separation in an annulus of physical radius r and width Ar 
at redshift zl- Moreover, contamination of the source sample 
by galaxies associated with the lens, or obscuration of back- 
ground galaxies by the optical counterpart of the lens, may 
lead to systematic error in the density measurement. A care- 
ful characterization of these errors is beyond the scope of this 
paper, whose purpose is to highlight that a simultaneous treat- 
ment of multiple lensing observables can lead to improved 
statistical performance and systematics self-calibration. 

However, some insight on the potential systematic impact 
of source clustering can be obtained by thinking of source 
clustering as a location-dependent modulation of the mean 
source density n. Consequently, removing the prior on h is 
roughly equivalent to the addition of source cluster noise. 
This scenario is shown in Figure |4] We see that there is an 
obvious degradation of the mass estimation, but even in this 
case a joint analysis is superior to a shear-only analysis. 

In light of this discussion, one might expect that the preci- 
sion that a real-world mass calibration experiment should be 



Albrecht, A., Bernstein. G., Cahn, R., Freedman, W. L., Hewitt, J., Hu, W., 
Huth, J., Kamionkowski, M., Kolb, E. W., Knox, L., Mather, J. C, Staggs, 
S., & Suntzeff, N. B. 2006, ArXiv Astrophysics e-prints, 
astro-ph/0609591 

Bartelmann, M. & Schneider, P. 2001, Phys. Rep., 340, 291 

Bernstein, G. M. 2009, ApJ, 695, 652 

Broadhurst, T. J., Taylor, A. N, & Peacock, J. A. 1995, Astrophys. J., 438, 
49 

Henry, J. P., Evrard, A. E., Hoekstra, H., Babul, A., & Mahdavi, A. 2009, 

ApJ, 691, 1307 
Jain, B. 2002, Astrophys. J., 580, L3 

Lombriser, L., Slosar, A., Seljak, U., & Hu, W. 2010, ArXiv: 1003.3009 
Mahdavi, A., Hoekstra, H., Babul, A., & Henry, J. P. 2008, MNRAS, 384, 
1567 

Mantz, A., Allen, S. W., Ebeling, H., & Rapetti, D. 2008, MNRAS, 387, 
1179 

MaiTone, D. P., Smith, G. P., Richard, J., Joy, M., Bonamente, M., Hasler, 
N., Hamilton-Morris, V., Kneib, J., Culverhouse, T., Carlstrom, J. E., 
Greer, C, Hawkins, D., Hennessy, R., Lamb, J. W., Leitch, E. M., Loh, 
M., Miller, A., Mroczkowski, T., Muchovej, S., Pryke, C, Sharp, M. K., 
& Woody, D. 2009, ApJ, 701, LI 14 

Meneghetti, M., Rasia, E., Merten, J., Bellagamba, F., Ettori, S., Mazzotta, 
P., Dolag, K., & Marri, S. 2010, A&A, 514, A93+ 



able to achieve lies somewhere between the dashed and dash- 
dotted line in Figure [4] Interestingly, one by-product of this 
kind of analysis is a high precision measurement of e. The 
posterior error in e from a joint analysis with known n, q, 
and a is nearly mass independent with a t w 0.04 assuming a 
source density n = 10 galaxies/arcmin 2 . This empirical deter- 
mination of e can be compared to its a priori expectation for 
a highly non-trivial test of residual PSF subtraction systemat- 
ics. Indeed, this is likely the most significant aspect of these 
type of analysis: we have explicitly demonstrated that by re- 
lying on multiple observables, it is possible to self-calibrate 
at least some subset of the parameters characterizing system- 
at ic uncertainties. A much more careful treatment akin to that 
of iBernsteinl (120091) is required to determine which systematic 
parameters are best calibrated using these methods. 

The degree of systematics-control is likely to be a deci- 
sive factor in upcoming surveys, as the statistical precision 
of shear-only mass calibration in upcoming photometric sur- 
veys is already exp ected to reach the 1 % level from standard 
shear-only analysis dRozo et al. 2010Q). Given that it remains 
to be seen whether systematics can be controlled well enough 
to saturate this limit, analyses which allow for some degree of 
self-calibration are doubly important: not only are the results 
automatically robust to the self-calibrated systematics, the 
empirical determination of systematic parameters can be com- 
pared to a priori expectations, allowing for crucial consistency 
tests of our understanding of the experiment. While we have 
not considered any systematics in the shear measurements 
here, we expect that the combination shear+density+size will 
be of comparable power in breaking degeneracies when in- 
cluding such sy s temat ics. Indeed, this is exactly the result of 
IVallinotto et al.l (120101) . These results provide strong motiva- 
tion to fully develop this type of joint lensing analysis as the 
main weak lensing mass calibration tool for near future data 
sets such as DES, Pan-Starrs, LSST, and Euclid. 
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